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The minimal model to describe many spin chain materials with ferroelectric properties is the 
Heisenberg model with ferromagnetic nearest neighbor coupling Ji and antiferromagnetic next- 
nearest neighbor coupling J2. Here we study the thermodynamics of this model using a density- 
matrix algorithm applied to transfer matrices. We find that the incommensurate spin-spin correla- 
tions - crucial for the ferroelectric properties and the analogue of the classical spiral pitch angle - 
depend not only on the ratio J2/I Ji| but also strongly on temperature. We study small easy-plane 
anisotropics which can stabilize a vector chiral order as well as the finite-temperature signatures 
of multipolar phases, stable at finite magnetic field. Furthermore, we fit the susceptibilities of 
LiCuV04, LiCu2C>2, and I^ZrCuCU. Contrary to the literature, we find that for LiCuVCU the best 
fit is obtained with J2 ~ 90 K and J2/I Ji| ~ 0.5 and show that these values are consistent with the 
observed spin incommensurability. Finally, we discuss our findings concerning the incommensurate 
spin-spin correlations and multipolar orders in relation to future experiments on these compounds. 

PACS numbers: 75.10.Jm, 75.40.Mg, 05.70.-a, 05.10.Cc 



I. INTRODUCTION 

A number of spin- 1/2 chain materials has recently been 
investigated in much detail which show multiferroic be- 
havior, i.e., an intricate interplay of magnetic and electric 
order. 1-7 A microscopic model for the electric ordering 
based on a spin current mechanism has been introduced 
in Ref. 8 and seems to provide an understanding for most 
of the experimental findings. Here the electric polariza- 
tion P is related to non-collinear spin-spin correlations, 
P <~ eij x (Si x Sj), where Si is a spin-1/2 operator at 
site i and a vector connecting the chain sites i and 
j. Later it has been shown based on a Ginzburg-Landau 
theory 9 and symmetry arguments 10 that this coupling 
between the polarization and magnetization always ex- 
ists independent of the crystal symmetry. To obtain 
a non-collinear spin structure on a lattice without ge- 
ometrical frustration, a minimal model has to contain 
additional interactions apart from the nearest neighbor 
Heisenberg interaction. One possibility are frustrating 
longer range interactions. For chain materials consisting 
of edge sharing copper-oxygen plaquettes like LICUVO4, 
LiCu202 or Li 2 ZrCuV04 the next-nearest neighbor in- 
teraction is particularly relevant leading to the minimal 
model 

H = ^ {Ji{S r S r+ \)xxz + J2{S r S r+ 2)xxz — hS z r } . 

r 

(i) 

Here (S r S r+l ) X xz = S?S? +i + Sy,S v r+l + AS z S z +l with 
A < 1 being an exchange anisotropy. For the edge shar- 
ing chains the nearest neighbor coupling is ferromagnetic 
( Ji < 0) while the next-nearest neighbor coupling is anti- 
ferromagnetic (J2 > 0). This is the case we want to study 
here. The external magnetic field is denoted by h. In the 
classical isotropic model (A — 1) without field the frus- 
tration leads to a helical spin arrangement with a pitch 



angle <j) = arccos(l/4a) for a = \ Jij J\ \ > 1/4. As in the 
classical model, the ground state of the quantum model 
(1) is ferromagnetic for a < 1/4. For a > 1/4 the ground 
state is a singlet whose nature has not fully been clari- 
fied yet. 11,12 Based on a renormalization group treatment 
starting from two decoupled Heisenberg chains (J\ = 0), 
it has been predicted that any small Ji < produces a 
finite but tiny excitation gap. 12 Numerically such a gap 
could not be resolved so far. 

Due to the SU(2) spin rotational symmetry only a 
quasi long-range helical order (algebraically decaying) is 
possible in the isotropic model (1) without field. Whether 
the correlation functions are indeed algebraically or in- 
stead exponentially decaying with a very large correla- 
tion length as suggested in Ref. 12 is a question which 
we will not address here and which is not important for 
the following discussions. One can in any case still de- 
fine a pitch angle <f> by studying the spin-spin correlation 
functions which, however, turns out to be substantially 
modified compared to the classical case due to quantum 
fluctuations. 13 If the SU(2) symmetry is broken by ap- 
plying either a magnetic field h or by an anisotropy A < 1 
then the vector chirality 

_ (g r x s r+n ) z = ^(SrS r+n — S r S^ +n ) (2) 

can have a nonzero expectation value because this re- 
quires only the breaking of the remaining Z 2 symmetry. 
The vector chirality as defined in (2) is directly related 

to the spin current j r = + 0^ which can be ob- 
tained from a continuity equation d t S z + d r j = 0. The 
equation of motion yields dS z = i[H,S z ] = — J\(k^ — 
~ ^(kt- 2 ^ — K i 2 \) which leads to the indentification 
j = J lK W+2J 2K ( 2 ). 14 

It has been shown recently that a magnetic field 
can also stabilize multipolar phases apart from a 
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phase 
well as 



with chiral order, both in the isotropic as 
in the anisotropic case. 14-17 Such multipo- 
lar phases are characterized by short-range transverse 
c x while correlations functions 



spin correlations (SqSj. 



(Sq S-^ • • • S r s, ; 



r+1 



are algebraically decaying in an 



?i-polar phase. 

In this article we want to investigate how the spin-spin 
correlations of the J1-J2 model (1) are affected not only 
by quantum but also by thermal fluctuations. At the 
end we want to relate our numerical results with vari- 
ous experiments on edge sharing copper-oxygen chains. 
A very powerful numerical method to study the thermo- 
dynamics of one-dimensional quantum chains directly in 
the thermodynamic limit is the density-matrix renormal- 
ization group applied to transfer matrices (TMRG). Here 
the Hilbert space is truncated to a small fixed number of 
states while the temperature is lowered successively. In 
the calculations presented here we will typically retain 
N = 60 — 200 states both in the system and the environ- 
ment block. We will present data for temperatures where 
the algorithm seems to be converged which is judged 
by comparing results obtained for different numbers of 
states. For a detailed description of the TMRG algo- 
rithm the reader is referred to Refs. 18-21. In a previous 
TMRG study the isotropic J\-J% model was investigated, 
however, correlation functions were not calculated. 22 

The manuscript is organized as follows: In Sec. II we 
study the isotropic model without magnetic field pay- 
ing special attention to the evolution of the incommen- 
surabilities with temperature. In Sec. Ill we consider 
the experimentally relevant case of a small easy-plane 
anisotropy. In Sec. IV we investigate signatures of mul- 
tipolar phases, which are stable for finite magnetic field, 
at finite temperatures. Finally, in Sec. V, we relate our 
numerical results to data for three experimentally well 
studied compounds (LICUVO4, LiCu202, Li 2 ZrCuV04). 
In particular, we fit the susceptibilities using model (1) 
which allows us to extract the parameters J\ and J 2 ■ We 
also discuss neutron scattering experiments which give 
access to the pitch angle and suggest future experiments 
regarding the realization and observation of multipolar 
orders. The last section is devoted to a brief summary 
and some conclusions. 



II. THERMODYNAMICS OF THE J x -J % CHAIN 

For the isotropic case, A = 1, without magnetic field 
it is known that there is a quantum critical point at 
a = 1/4 separating a ferromagnetic phase for a < 1/4 
from a phase with a singlet ground state. A tiny gap for 
a > 1/4 has been predicted, 12 however, even if this gap 
really exists it will only show up at extremely low temper- 
atures. In the temperature range accessible numerically, 
the susceptibility x f° r a ~ > 00 smoothly approaches the 
result for the nearest-neighbor Heisenberg chain known 
from Bethe ansatz 23 as is shown in Fig. 1. At the crit- 
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FIG. 1: (color online) (a) Susceptibilities near the quan- 
tum critical point a — 1/4. (b) Susceptibilities for a = 
2.0, 1.6, 1.2, 1.0, 0.8, 0.6 (solid lines from bottom to top). For 
comparison, the result for the nearest-neighbor Heisenberg 
chain obtained by Bethe ansatz is shown (dashed line). 




FIG. 2: (color online) Specific heat for a = 0.4,0.6, •■• ,2.0 
(from bottom to top). 



ical point, x(T) will show a power-law divergence. 24 A 
thorough investigation of the critical properties using the 
TMRG algorithm will be presented elsewhere. 25 

Numerically, we can also obtain the inner energy e and 
the free energy / (see Ref. 21). This gives us the entropy 
s — (e— f)/T and the specific heat C = Tds/dT using a 
numerical derivative. The results for various a are pre- 
sented in Fig. 2. For all frustration parameters we find 
a broad maximum which shifts monotonically to larger 
temperatures T / J2 with increasing a. 

Next, we want to study various correlation functions. 
At finite temperatures all correlation functions will be 
exponentially decaying and we can expand any two-point 
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FIG. 3: (color online) The wave vector k of the spin-spin 
correlation function depends not only on a but also on tem- 
perature. With increasing temperature, k first increases for 
a close to a c while it decreases for larger a. 

correlation function of an operator O r as 

oo 

(O O r ) - (Oo) (O r ) = Mje- r ^ cosfor) . (3) 

3=0 

Here Mj is a matrix element, £j a correlation length, and 
kj the corresponding wave vector. Within the TMRG 
algorithm both £j and kj are determined by the ratio 
of the leading to next-leading eigenvalues of the transfer 
matrix. 18,21 For large distances r > 1 the correlation 
function is dominated by the largest correlation length 
and the corresponding wave vector which we will denote 
as £ = £o and k = kg in the following. 

Of particular interest is the evolution of the wave vec- 
tor k for the spin-spin correlation (SoS r ). This can be 
seen as the quantum mechanical analogue of the pitch 
angle of the classical spiral. In Fig. 3 it is shown that 
k does not only depend on the frustration a but also 
shows a strong dependence on temperature, in particu- 
lar, for values of a close to the critical point a c = 1/4. 
Contrary to the classical case, the pitch angle at low tem- 
peratures is very close to ir/2 for a > 0.6 as is shown 
in more detail in Fig. 4(a). This means that for these 
frustration values the spin structure develops an almost 
perfect (short-range) four-site periodicity at low temper- 
atures. In Fig. 4(b) the extrapolated k value (pitch an- 
gle) for zero temperature is shown as a function of a. 
Our result is in good agreement with an earlier DMRG 
calculation. 13 We also find that with increasing temper- 
ature the pitch angle becomes again more classical for a 
frustration a > 0.5 (see squares in Fig. 4(b)). 

One might think that the strong trend towards a for- 
mation of a four-site periodic structure with increasing 
a could be related to a concomitant dimerization. How- 
ever, as is shown exemplarily in Fig. 5, this is not the 




FIG. 4: (color online) (a) k as a function of temperature for 
a = 0.6, 0.8, • • • , 2.0 (from bottom to top), (b) Extrapolated 
k (pitch angle) for zero temperature (circles), and k at T / J2 = 
1.0 (squares). The dashed lines are a guide to the eye. For 
comparison, the classical pitch angle (j> = arccos(l/4a) (solid 
line) is shown. 
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FIG. 5: (color online) Various correlation functions for a = 2 
and T/J2 = 0.05. The circles denote the spin-spin (O r = S r ), 
the squares the dimer (O r = S r S r +i), and the diamonds the 
chiral correlation function (O r — S r X S r +i)- The lines are a 
guide to the eye. 



case. While substantial dimer and chiral correlations 
with comparable correlation lengths do exist, both cor- 
relation lengths are about a factor 2 smaller than the 
spin-spin correlation length. 

Another way of defining a quantum analogue of the 
classical pitch angle is the possibility to study at which 
wave vector q the static spin structure factor 



S(q) 



2^2cos(qr)(S S r 



(4) 



r=l 



is peaked. This definition of the pitch angle is equiva- 
lent to the definition in (3) if the correlation function is 
dominated by the leading correlation length, i.e., £0 3* 
for j > 1. This is the case for all frustration parameters 
shown in Fig. 6 except for the critical value a c = 1/4 
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FIG. 6: (color online) The static spin structure factor S(q) for 
temperatures T/ 1 Ji | = 0.5,0.2,0.1,0.05,0.035,0.02 (along the 
arrow direction) with (a) a = 1/4 (quantum critical point), 
(b) a = 0.3, (c) a = 0.6, and (d) a = 1.0. The inset of (a) 
shows the deviation from the sum rule S = ~}2 S(q) = 3/4. 



(Fig. 6(a)). Here the structure factor is very fiat near 
q ~ meaning that in the expansion (3) comparable cor- 
relation lengths exist with a wave vector k = and with 
wave vectors which have small incommensurate values. 
The accuracy of the numerical data can be checked by 
calculating the sum rule. For all a and all temperatures 
the sum rule is fulfilled with deviations of the order 1CP 3 
to 1CP 4 as is exemplarily shown in the inset of Fig. 6(a). 



III. THE ANISOTROPIC CASE 

In the edge-sharing spin chain compounds substan- 
tial exchange anisotropics exist. For LiCuVC>4, for ex- 
ample, the g tensor has been measured by ESR and 
an exchange anisotropy of the order of 6% percent has 
been estimated. 26 This anisotropy is expected to pin the 
spins to the ab plane and this is indeed observed in 
experiment. 5,7 Theoretically, one expects that an XX Z- 
type anisotropy (A < 1) enhances chiral correlations. In 
this case even long range chiral order ((«("') ^ 0) is pos- 
sible at zero temperature in the purely one-dimensional 
model because only the remaining Z2 symmetry has to 
be broken. In previous numerical studies the anisotropic 
model (1) has already been investigated, the results, how- 
ever, have been contradictory 17,27 In Ref. 27 a large re- 
gion of the phase diagram has been found to be occu- 
pied by a dimer phase whereas spin liquid phases oc- 
cur for very small and large values of a. In Ref. 17 
a dimer phase and spin liquid phases have again been 
found but in addition also a chiral phase for a ^> 
1. Both studies were based on exact diagonalization. 
In a very recent density-matrix renormalization group 
study the anisotropic model at finite magnetization was 
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FIG. 7: (color online) The longitudinal spin (solid line), dimer 
(dashed line), and chiral correlation length (dot-dashed line) 
for a — 0.4 and A — 0.8. At low temperatures the chi- 
ral correlation length diverges stronger than 1/T indicating 
long-range order at zero temperature, (a) The longitudinal 
spin-spin and the chiral correlation function at T/J2 = 0.125 
showing that chiral correlations dominate, (b) The wave vec- 
tor k of the longitudinal spin-spin correlation function. 

investigated. 15 A dimer phase was not found to be stable, 
instead it was shown that for most parameters a chiral 
phase or a spin liquid phase (SDW2 phase) occur. 

We will first concentrate on one value of anisotropy, 
A = 0.8. This is larger than what is expected for the 
edge- sharing cuprate chains. However, this way the ef- 
fects of the anisotropy are more obvious and the results 
should qualitatively be similar to those for realistic val- 
ues for these materials. We will study the case of zero 
field to see whether a dimer or a chiral phase occur. In 
Fig. 7 results for various correlation lengths at a = 0.4 
are shown. For correlation functions which decay alge- 
braically at zero temperature we expect that the correla- 
tion length diverges as £ ~ 1 jT. Studying £T as in Fig. 7 
thus separates long and short ranged correlations. In this 
case we see that the chiral correlation length diverges 
stronger than 1/T indicating that for these parameters 
we do have long range chiral order at zero temperature. 
While the dimer is indeed larger than the chiral correla- 
tion length over a wide temperature range, the situation 
is reversed at low temperatures. Hence a dimer phase 
can only possibly occur if also interchain couplings are 
present which might stabilize such an order at interme- 
diate temperatures. That the chiral correlations indeed 
dominate at low temperatures is shown exemplarily for 
T/J2 = 0.125 in Fig. 7(a) where the chiral and longitu- 
dinal spin-spin correlation functions are compared. Here 
the wave vector k of the longitudinal spin-spin correlation 
function is incommensurate and again shows a strong de- 
pendence on temperature (Fig. 7(b)). 

For a = 0.6, shown in Fig. 8(a), the chiral correla- 
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FIG. 8: (color online) The longitudinal spin (solid line), dimer 
(dashed line) and chiral (dot-dashed line) correlation lengths 
for A = 0.8 with (a) a = 0.6 and (b) a = 2.0. 



tion length is again dominant at the lowest temperatures 
which are accessible numerically, however, there is no in- 
dication for long range chiral order at zero temperature. 
Remarkably, all three correlation lengths are of very sim- 
ilar magnitude at low temperatures. This suggest that 
we are very close or at the phase transition from a phase 
with long range chiral order at smaller values of a to a 
phase which most likely has spin liquid character with 
algebraically decaying correlation functions. Indeed, at 
a = 2.0 shown in Fig. 8(b) the spin-spin correlation 
length is clearly dominant with oscillations k ss ir/2 at 
low temperatures pointing to an SDW2 phase. We con- 
clude that if a dimer phase exist at all for A = 0.8 it 
has to be confined to a very narrow range of frustration 
parameters 0.4 < a < 0.6. In Ref. 27 an unidentified 
phase was found for 0.3 < a < 0.45 with symmetries as 
expected for the chiral phase. For a > 0.7, on the other 
hand, a spin liquid phase was predicted with a dimer 
phase in between these two phases. In our TMRG calcu- 
lations we never find dominant dimer correlations at low 
temperatures and our results suggest that a dimer phase 
might not be stable at all but instead a direct phase tran- 
sition from the chiral to a SDW 2 phase atasj 0.6 occurs. 



IV. FINITE MAGNETIC FIELD: MULTIPOLAR 
PHASES 

In the multiferroic spin chain compounds a magnetic 
field can be used to switch the electric polarization. A 
sufficiently strong field induces a flop of the spins from the 
easy-plane spiral to a spiral perpendicular to the applied 
magnetic field. According to the spin current mechanism 
this also leads to a rotation of the electric polarization. 
Such a switching of the ferroelectric polarization by an 
applied magnetic field has been observed in LiCuV04 7 
as well as in LiCu202- 2 Numerical studies of the J1-J2 




FIG. 9: (color online) The longitudinal spin (solid line), dimer 
(dashed line) and chiral (dot-dashed line) correlation lengths 
for a = 0.4 with A = 1 and h = 0.25 J2. (a) The wave vector 
k of the spin-spin correlation, (b) The magnetization m as a 
function of temperature. 

model have shown that small magnetic fields stabilize 
the chiral order while multipolar phases become stable at 
intermediate field strengths before the system ultimately 
becomes fully polarized for large fields. 14,16,17 

As already eluded to in the introduction, an algebraic 
decay of the (SqS^ ■ ■ ■ S~ S~ +1 • • •) correlation function is 



expected in an n-polar (n > 2) phase while the transverse 
spin correlation function is gapped. For the oscillations 
of the longitudinal spin correlation function we expect in 
general 

fcr^o = -(l-2m) (5) 
n 

where m is the magnetization and is directly related to 
the shift of the Fermi points. For zero field, n — 2 would 
correspond to the 4-site periodic spin structure discussed 
previously. 

An particular interesting case is a — 0.4 with h = 
0.25 J2 shown in Fig. 9. Here we find that the longitu- 
dinal spin correlation length is largest at high tempera- 
tures, then there is a temperature range where the chiral 
correlation length dominates while at the lowest acces- 
sible temperatures the spin correlation again dominates 
and seems to diverge like 1/T. For zero temperature this 
field corresponds to a magnetization m 0.15 ~ 0.3m sa t 
(Fig. 9(b)) where m sa t = 0.5 is the saturation magne- 
tization. Comparing with the zero temperature phase 
diagrams obtained in Refs. 14,16 we see that for these 
values we are in the SDW2 phase but very close to the 
phase boundary with the chiral phase. Our results ba- 
sically seem to confirm this picture although the oscil- 
lations k apparently slightly deviate from the value ex- 
pected for a n = 2 multipolar phase even at the lowest 
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FIG. 10: (color online) S zz (q) for a = 0.32, h = 0.02 Ji and 
temperatures T/\.h\ = 2.0,1.0,0.5,0.2,0.1,0.07,0.05 (from 
botom to top). The dashed line marks q/[ir(l — 2m)] = 1/3. 
The inset shows the magnetization m as function of temper- 
ature. 



temperatures (see Fig. 9(a)). Our data also show that 
in a certain temperature window chiral correlations can 
still be dominant. This might be relevant once inter- 
chain couplings are taken into account and might lead to 
a chiral phase stable at intermediate temperatures while 
the spins are collinear at higher and lower temperatures. 
Two magnetic phase transitions have indeed been ob- 
served in LiCu202 where first a sinusoidal magnetic or- 
der is established followed by a helical order at lower 
temperatures. 2,28 While the magnetic structure even in 
the low temperature phase apparently is much more com- 
plicated than a simple planar spiral 28 and the J1-J2 
model does not seem to capture the essential physics of 
this compound (see next section) it is nevertheless inter- 
esting that even in this very simple model phases might 
be stable only in a certain temperature range thus pos- 
sibly giving rise to multiple magnetic phase transitions 
once interchain couplings are taken into account. 

In Fig. 10 the longitudinal spin structure factor S zz (q) 
is shown for a — 0.32 and h = 0.02Ji. At low tempera- 
tures S zz (q) is peaked at q = tt(1 — 2m) /3 (dashed line 
in Fig. 10) indicating an n = 3 multipolar phase. Thus 
experimentally a multipolar phase can already be inden- 
tified at finite temperatures by studying the evolution 
of the structure factor as obtained in neutron scattering. 
Such an experiment would be most interesting for a com- 
pound close to the quantum critical point a c — 1/4 where 
small magnetic fields are sufficient to stabilize n — 3 or 
even n = 4 multipolar phases. 16 As discussed in detail in 
Ref. 4 and in the following section, Li 2 ZrCu04 seems to 
be a very promising candidate for such a material. 



V. MULTIFERROIC SPIN CHAIN MATERIALS 

While in the corner sharing copper-oxygen chain com- 
pounds like Si^CuOa and SrCu02 the antiferromagnetic 
exchange interaction is of very similar magnitude, J ~ 
2000 K, a wide range of exchange parameters J\ and 
J2 can be found in the literature for the edge sharing 
compounds. For LiCuVC>4, for example, fits of the sus- 
ceptibility and of neutron scattering data have led to the 
estimate J\ ~ —20 K and J2 ~ 50 K so that a ~ 2. 5. 5 
For Li2ZrCu04, on the other hand, susceptibility and 
specific heat data have been fitted by using J x ~ —300 
K and a ~ 0.3. 4 Given that the chains in these com- 
pounds consist of the same edge sharing copper-oxygen 
plaquettes, this huge variation in the magnitude of the 
exchange couplings is surprising. 

Here we want to reanalyze the data of three of 
the best studied multiferroic chain compounds namely 
Li 2 ZrCu04, LiCuVC>4 and LA2CUO2. We will concen- 
trate on fitting susceptibility data using 

Xoxp = X0+XJ1-J2 ■ ( 6 ) 

Here Xcxp is the experimentally measured susceptibil- 
ity Xo is a constant contribution due to core diamag- 
netism and Van-Vleck paramagnetism and XJ1-J2 is the 
numerically calculated susceptibility for the J\-Ji model 
(1). Such fits work extremely well for the corner shar- 
ing compounds because the intrachain coupling is about 
three orders of magnitude smaller than the interchain 
coupling. 29-33 For LiCuV04 it has been reported that a 
three-dimensional (3D) magnetic order becomes estab- 
lished below Tj,b ~ 2.3 K. 6 While this points to intra- 
chain couplings which are only one or at most two orders 
of magnitude smaller than the interchain couplings, a 
purely one-dimensional model is still expected to be a 
good approximation as long as T 3> T$d. In L^ZrCuO^ 
the situation has not fully been clarified yet. A possible 
phase transition might occur at T ~ 6 K. 4 Even if this is a 
3D ordering transition, a one-dimensional model should 
again be valid over a wide temperature range. There- 
fore the expectation is that the physics of LiCuV04 and 
Li 2 ZrCu04 can largely be understood within the frame- 
work of the 3\-Ji model with the spin-current mechanism 
being responsible for the multiferroic properties. 

The situation for LiCu 2 02, on the other hand, seems to 
be much more involved. 1 ~ 3 ' 28,34 Despite a number of stud- 
ies, the magnetic structure remains controversial. For 
instance, both a spiral spin order in the ab as well as 
in the be plane have been reported. 1 ' 3 If the spin cur- 
rent model is the correct explanation for the observed 
polarization along the c-axis 2 then the be spiral must be 
realized. The two magnetic ordering transitions at com- 
paritatively large temperatures, T 22 K and T as 24 
K, 28 point to much larger interchain couplings than for 
the two other compounds discussed above. In fact, it has 
been found that there are substantial spin correlations 
not only along the chain direction (5-axis) but also in 
the plane of the copper-oxygen plaquettes perpendicular 
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to the chain direction (a-axis) making the compound at 
low temperatures almost two dimensional. 34 One might 
speculate that the reason why the magnetic properties of 
this material are so different is related to the fact that 
two different copper sites exist - the in-chain Cu 2+ and 
the Cu 1+ interconnecting different chains. This might 
lead to substantial charge fluctuations and thus to en- 
hanced magnetic interchain couplings. Nevertheless, a 
fit of the susceptibility at high temperatures using the 
J\-Ji model might still be useful to obtain an estimate 
for the magnitude of the J\ and J 2 couplings. 

The absolute values of the measured maxima of x al- 
ready allow some general statements about the magni- 
tude of the exchange constants and the frustration pa- 
rameters a. In LiCuV04 and in LiCu202 the suscepti- 
bility is about two orders of magnitude larger than in the 
corner sharing chain compounds Si^CuOa and SrCuC>2 so 
that the antiferromagnetic exchange constant J2 should 
also roughly be two orders of magnitude smaller. In 
Li2ZrCu04 the maximum is almost an order larger than 
in the other two compounds clearly indicating that this 
compound should be closer to the quantum critical point 
a c = 1/4 than the other two compounds. If the mag- 
netic exchange constants are of comparable magnitude, 
then a is expected to be largest for LiCu202 and small- 
est for Li 2 ZrCu04 with LiCuV04 having an intermediate 
frustration parameter. 

The susceptibility of LiCuVC>4 shown in Fig. 11 has 
a maximum and a local minimum at low temperatures. 
If x can indeed be described by the Jy J 2 model then a 
comparison with Fig. 1 indicates that a ~ 0.5. The large 
a used in Ref. 5 can certainly not explain this structure. 
By performing a fit according to (6) we find that J2 = 91 
K and a = 0.5 yields the best fit as shown in Fig. 1(a). 
A good fit is also possible with J2 — 72 K and a = 
0.6 (Fig. 1(b)). If we reduce the next-nearest neighbor 
exchange to J2 ~ 50 K as assumed in Ref. 5 we have to 
choose a ~ 1 to obtain the best fit possible with this value 
of J2 (Fig. 1(c)). However, such a fit fails to reproduce 
the low temperature structure. 

The susceptibility data for Li 2 ZrCu04 have already 
been analyzed in Ref. 4 by comparing with exact diago- 
nalization data for small rings consisting of up to N — 20 
sites. Finite size effects are expected to be neglegible if 
T 3> v/N where v is the spin velocity. The spin velocity 
is of the order of the exchange constant so that the finite 
size data for x should be reliable for T > 20 K. Given 
that the maximum of x is located at T « 7.6 K, i.e., at 
temperatures where finite size effects might play a role, 
it is helpful to reanalyze the experimental data using the 
TMRG algorithm which yields results directly for the in- 
finite system. As shown in Fig. 12 very similar fits to the 
one in Ref. 4 are obtained. While the magnitude of the 
frustration parameter can only be varied slightly if one 
wants to obtain a reasonable fit, the exchange parame- 
ters change dramatically, for example, from J\ = —273 
K for a = 0.3 to J x = -75 K for a = 0.4. In LDA calcu- 
lations Ji — —151 ± 35 K and J2 = 35 ± 12 K was found 
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FIG. 11: (color online) The experimentally measured sus- 
ceptibility of LiCuVC>4 for H \\ c (solid line) compared 
to fits (dots) following Eq. (6) with a gyromagnetic ratio 
g = 2.313. 6 ' 26 (a) The best fit is obtained with Xo = 6 ■ 1(T 5 
emu/mol, J2 = 91 K, and a — 0.5. (b) An alternative fit with 
Xo = 1.4 ■ 10" 4 emu/mol, J 2 = 72 K, and a = 0.6. (c) An ex- 
change coupling J2 = 52 K similar to the one in Ref. 5 can be 
used with xo = 2.7 • 10 -4 emu/mol, however, the frustration 
parameter a — 1.0 chosen to obtain the best fit is still much 
smaller than the one in Ref. 5 and the fit is much worse than 
the ones shown in (a) and (b). 
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FIG. 12: (color online) The measured susceptibility (solid 
line) for Li2ZrCu04 taken from Ref. 4 compared to fits using 
the Ji — J2 model. An excellent fit (circles) is obtained using 
a = 0.3, g = 2 and Ji = -273 K (J 2 = 81.9 K) with xo = 
confirming the analysis in Ref. 4 based on exact diagonal- 
ization data. A reasonable fit (squares) is also possible with 
a = 0.4, g = 2.2, xo = 0, and Ji = —75 K. Inset: Blow-up of 
the low temperature region. 



which also will allow for a reasonable fit with a frustra- 
tion parameter a € [0.3, 0.4]. 4 We therefore confirm the 
main conclusion of Ref. 4 that Li2ZrCu04 is close to the 
critical point a c = 1/4. 

Finally, we want to analyze the susceptibility data for 
LiCu202- As already mentioned, we do not expect that 
the J1-J2 model will describe the susceptibility of this 
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FIG. 13: (color online) The measured susceptibility (solid 
line) for LiCu2 02 with H \\ c taken from Ref. 1 compared to 
fits using the Ji — J2 model. With g — 2.2 the best fit is 
obtained with a — 0.6, J2 = 120 K and xo = (circles) while 
for g — 2.3 a fit with a = 1.6, J2 = 83 K and xo = works 
best (squares). 

compound at low temperatures and indeed we find that 
it is impossible to obtain a good fit down to low tem- 
peratures (see Fig. 13). If we concentrate on temper- 
atures large compared to the magnetic ordering transi- 
tions at T ~ 20 K, then a fit is possible and we find that 
J 2 ~ 80 — 120 K. Because the susceptibility can only be 
fitted at high temperatures, the frustration parameter 
cannot be fixed precisely. 

It is, however, important to note that the best fits for 
all three compounds yield values of J2 < 100 K while 
J l r*, —200 K. This confirms our expectation that these 
materials consisting of the same edge sharing copper- 
oxygen plaquettes should have very similar magnetic ex- 
change constants as is also the case for the corner shar- 
ing chain compounds. The values we obtain from the 
susceptibility fits are consistent with values in the liter- 
ature for LiCu 2 2 (Ji 120 K, J 2 - 80 K) 1 ' 2 and 

for Li 2 ZrCu04. 4 For LiCuV04, however, we find values 
which differ dramatically from the values given in Ref. 5 
which later were also used in a number of other publica- 
tions. We want to stress again that these values are not 
consistent with the latest susceptibility data. 6 An anal- 
ysis of neutron scattering data using a standard spin- 
wave dispersion which seems to confirm these values 5 is 
in our opinion not applicable here. In such an analysis 
the magnon bandwidth is directly determined by the bare 
exchange couplings while the stark deviation of the quan- 
tum from the classical pitch angle (see Fig. 4(b)) suggests 
that the frustration parameter is strongly renormalized 
due to quantum fluctuations. 

The pitch angle measured experimentally in Ref. 5 of 
cj) « 83.7° is in fact in excellent agreement with the frus- 
tration parameters a — 0.5 — 0.6 obtained from the fits 
in Fig. 11. According to Fig. 4(b) we have a quantum 
pitch angle (j> ss 83.5° for a = 0.5 and cf> w 87.3° for 



a = 0.6. It is important to stress again that only in the 
classical model large frustration parameters are needed 
to obtain pitch angles close to 90°. We expect, accord- 
ing to Fig. 4(a), that the pitch angle in LiCuVC>4 can 
be reduced by 15 — 20% by increasing the temperature 
to T ~ J 2 ~ 90 K. For Li 2 ZrCuC>4 the magnetic struc- 
ture has not been studied so far. The pitch angle for 
a = 0.3 at zero temperature obtained from our numerical 
calculations is <fi « 51.5°. Here we expect a large varia- 
tion with temperature (see Fig. 3 and Fig. 6(b)) and it 
would be interesting to see if this can also be observed 
experimentally. For LiCu 2 2 a pitch angle <j> « 62.6° 
has been measured. 1 ' 2 Such a small pitch angle cannot 
be explained within the J\- J 2 model given that a > 0.6 
according to the fits shown in Fig. 13. We are therefore 
again lead to the conclusion that the J1-J2 model cannot 
explain the experimental data for this compound. 



VI. CONCLUSIONS 

The rich physics of the J\-J% model with ferromagnetic 
nearest neighbor coupling J\ and antiferromagnetic next- 
nearest neighbor coupling J 2 has attracted a lot of inter- 
est recently. The phase diagram for this model includ- 
ing exchange anisotropies and magnetic fields has been 
addressed in a number of numerical studies. 13-17 Here 
we have shown that the physical properties of this sim- 
ple model are even more intriguing if the interplay of 
quantum and thermal fluctuations is taken into account. 
In particular, we found that the incommensurate oscil- 
lations of the spin-spin correlation function (the quan- 
tum analogue of the pitch angle of the classical spiral 
order) does not only strongly depend on the frustration 
parameter a — J 2 /|Ji| but also on temperature. For 
zero temperature we find an incommensurability (pitch 
angle) <\> f° r a ^ 0-6 which is very close to 90° in the 
quantum model and thus very different from the classical 
pitch angle (f> — arccos(l/4a) in agreement with an ear- 
lier study. 13 At temperatures T ~ J 2 , however, the pitch 
angle is much closer to the classical value for a > 0.5. 
Furthermore, we find a very strong temperature depen- 
dence of the pitch angle for frustration parameters close 
to the critical point a c = 1/4. We therefore expect that 
the wave vector where the static spin structure factor for 
Li 2 ZrCu04 is peaked - a compound which according to 
Ref. 4 and the susceptibility analysis performed here has 
a frustration parameter a ~ 0.3 — 0.4 - varies significantly 
with temperature. For this range of frustration parame- 
ters we find that a small easy-plane anisotropy (at zero 
magnetic field) leads to long-range chiral order even in 
the purely one-dimensional model. 35 We could, however, 
find no evidence for the dimcr phase which was predicted 
in Ref. 27 for the anisotropic case for larger frustration 
parameters. Another observation which again might be 
relevant for future studies on Li 2 ZrCuC>4 is that small 
magnetic fields (which would correspond to 5 — 7 T for 
Li 2 ZrCuC>4) can stabilize an SDW3 (n = 3 multipolar) 
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phase for a ~ 0.3. We suggest that such a phase can 
already be identified by neutron scattering at finite tem- 
peratures by monitoring at the same time the peak posi- 
tion g max of the static longitudinal structure factor and 
the magnetization m. The signature of the SDW3 phase 
would be that <7 max — > ir(l — 2m) /i for T — > 0. 

We also reanalyzed the susceptibility data for 
LiCuVC>4. Here the newer data in Ref. 6 seem to be of 
much better quality than the older data in Ref. 5. The 
newer data clearly show that the exchange parameters 

Ji 20 K, J 2 - 50 K and a = 2.5 found in Ref. 5 

do not yield a reasonable fit. These exchange parameters 
also seem very unlikely given that they deviate signif- 
icantly from those in other edge-sharing copper-oxygen 
chains. Our analysis shows that the susceptibility data 
are best fitted with J 2 w 70 - 90 K and a = 0.5 - 0.6. 

In LiCu202 two magnetic ordering transitions already 
occur at temperatures T^d ~ 20 K and an analysis based 
on the purely one-dimensional Ji-J 2 model clearly can- 
not be as successful as for LiCuVC>4 and Li2ZrCu04. 
For T 3> Tsd we have shown that a reasonable fit 
of the susceptibility is nevertheless possible leading to 
J 2 ~ 80 — 120 K. The frustration parameter remains, 
however, ambiguous in such a high temperature fit and 



we find a € [0.6, 1.6] with the best fit being obtained for 
a = 0.6. 

Based on the best fits of the susceptibilities we con- 
clude that all three considered compounds seem to have 
very similar exchange constants J\ ~ —200 K and J 2 < 
100 K. The frustration parameter, on the other hand, 
varies from a ss 0.3 for Li 2 ZrCu04, a = 0.5 — 0.6 for 
LiCuV04 toa> 0.6 for LiCu 2 02- The parameters found 
here for LiCuV04 are fully consistent with the measured 
pitch angle ~ 84°. For LiCu 2 C>2, on the other hand, 
the small measured pitch angle ~ 63° cannot be ex- 
plained within the J1-J2 model stressing again that this 
model is not sufficient to explain the experimental data 
for this compound. The smallest pitch angle is expected 
for Li 2 ZrCu04. Based on our numerical calculations we 
predict <fi <~ 53°. 
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